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Next-to-leading order corrections to jet cross sections in deep inelastic scattering 
at HERA are studied. The predicted jet rates allow for a precise determination of 
a s (fJ.n) and the gluon density at HERA. We argue, that the "natural" renormal- 
ization and factorization scale is in general set by the average of the jets in the 
Breit frame. Some implications for the associated forward jet production in the 
low x regime at HERA are discussed. 



1 Introduction 

One of the main topics to be studied at HERA is the deep inelastic production 
of multi jet events- where good event statistics are expected allowing for preci- 
sion tests of QCDMeI. Clearly, next-to-leading order (NLO) QCD corrections 
are mandatory on the theoretical side for such tests. Full NLO corrections for 
one and two-jet production cross sections and distributions are now available 
□ and implemented in the fully differential ep — > n jets event generator MEP- 
JET, which allows to analyze arbitrary jet definition schemes and general cuts 
in terms of parton 4-momenta. A variety of topics can be studied with these 
tools. They include: a) The determination of a s (fj,n) over a range of scales \xr 
from dijet production, b) The measurement of the gluon density in the proton 
(via 717 — > qq) c) The measurement of the polarized gluon density Ag with 
polarized beams of electrons and protons at HERA El. d) Associated forward 
jet production in the low x regime as a signal of BFKL dynamics. 

The importance of higher order corrections and recombination scheme de- 
pendencies of the two jet cross sections for four different jet algorithms (cone, 
kx, JADE, W) was already discussed in Refs. [4,6,7]. While the higher or- 
der corrections and recombination scheme dependencies in the cone and kx 
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schemes are small, very large corrections can appear in the W-scheme. De- 
pending on the definition of the jet resolution mass and on the recombination 
schem£,_the NLO cross sections in the VF-scheme can differ by almost a factor 
of twoLTtl TrefzgerEl and RosenbauerB find similarly large differences in the ex- 
perimental jet cross sections (which are in good agreement with the theoretical 
predictions from MEPJET), if the data are processed with exactly the same 
jet resolution mass and recombination prescription as used in the theoretical 
calculation. Note however, that the dependence on the recombination scheme 
first appears in the NLO calculation, where a jet may be composed of two par- 
tons. This internal jet structure is thus only simulated at tree level and thus 
the dependence of the cross section on the recombination scheme is subject 
to potentially large higher order corrections. An alternative fully differential 
NLO MonteCarlo program DISENT is under construction 

Previous programstj£3 were limited to a JADE type algorithm and are not 
flexible enough to implement the various recombination schemes. In addition, 
approximations were made to the matrix elements in these programs which 
are not valid in large regions of phase space □. It is therefore nat surprising, 
that one finds inconsistent values for a s using these programs □, if data are 
processed with different recombination schemes. 

We conclude from these studies that the cone and kr schemes appear better 
suited for precision QCD tests and will concentrate only on those schemes in 
the following. 

In the cone algorithm (which is defined in the laboratory frame) the dis- 
tance AR = [Arf) 2 -I- (A0) 2 between two partons decides whether they 
should be recombined to a single jet. Here the variables are the pseudo-rapidity 
r\ and the azimuthal angle <f>. We recombine partons with AR < 1. Further- 
more, a cut on the jet transverse momenta of pr(J) 5 GeV in the lab frame 
is imposed. 

For the kr algorithm (which is implemented in the Breit frame) , we follow 
the description introduced in Ref. [11]. The hard scattering scale Ej, is fixed to 
40 GeV 2 and y cut — 1 is the resolution parameter for resolving the macro-jets. 



2 The determination of a s (nn) from dijet production 

The dijet cross section is proportional to a s (/i^) at leading order (LO), thus 
suggesting a direct measurement of the strong coupling constant. However, 
the LO calculation leaves the renormalization scale (1r undetermined. The 
NLO corrections substantially reduce the renormalization and factorization 
scale dependencies which are present in the LO calculations (see below) and 
thus reliable cross section predictions in terms of a s (mz) are made possible. 
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Clearly, a careful study of the choice of scale in the dijet cross section is needed 
in order to extract a reliable value for a s (M z ). Jet production in DIS is a 
multi-scale problem and_.it is not a priori clear at which scale a s is probed. 
However, it was argued 0, that the "natural" scale for jet production in DIS 
is set by the average k^ of the jets in the Breit frame. Here, (k^(j)) 2 is 
defined by 2Ej(l — cosOjp), where the subscripts j and P denote the jet and 
proton, respectively (all quantities are defined in the Breit frame). It can 
be shownLJ that £V U) smoothly interpolates between the correct limiting 
scale choices, it approaches Q in the parton limit and it corresponds to the jet 
transverse momentum p 8 , (with respect to the 7*-proton direction) when the 
photon virtuality becomes negligible. It therefore appears to be the "natural" 
scale for multi jet production in DIS. 

Fig. la shows the scale dependence of the dijet cross section in LO and 
NLO for the kr scheme. We have considered scales related to the scalar sum 
of the partons (solid curves), the scalar sum of the partons pr with respect 
to the 7*-boson direction (dashed curves) and the virtuality Q of the incident 
photon. The LQ-/NLO) results are based on the LO (NLO) parton distribu- 
tions from GRVE-- together with the one-loop (two-loop) formula for the strong 
coupling constant. Kinematical cuts are imposed to closely model the HI event 
selection. More specifically, we require an energy cut of E(e') > 10 GeV on the 
scattered electron, and a cut on the pseudo-rapidity 77 = — lntan(0/2) of the 
scattered lepton depending on Q 2 , i.e. -2.794 < r](l') < -1.735 for Q 2 < 100 
GeV 2 -1.317 < r](l') < 2.436 for Q 2 > 100 GeV 2 . In addition, we require 
— 1.154 < r)(j) < 2.436. The scale dependence of the dijet cross section does not 
markedly improve in NLO for fi 2 = £Q 2 (dotted lines in Fig. la). The resulting 
£ dependence for = n% = £ (YV k^(i)) 2 (fi 2 R = /i| = £ (£. p* (i)) 2 ) is 
shown as the solid (dashed) lines in Fig. la. The uncertainty from the varia- 
tion of the renormalization and factorization scale scale is markedly reduced 
compared to the LO predictions. Hence, the theoretical uncertainties due to 
the scale variation are small suggesting a precise determination of a s (< k^ >) 
for different < k§ > bins, where < kj, >— | {Y2j=i 2 ^tU))- ^ has also been 
shown 0, that the < kr >, < Pt > and Q distributions can differ substan- 
tially for different < k^ > bins (in particular for low Q 2 ), which explains the 
large differences in the scale dependence in Fig. 1. Fig. lb shows the scale 
dependence of the dijet cross section in cone scheme for cuts similar to the 
ZEUS selection: 40 GeV 2 < Q 2 < 2500 GeV 2 , 0.04 < y < 1, and a cut on the 
pseudo-rapidity 77 = — lntan(0/2) of the scattered lepton and jets of \rj\ < 3.5. 
The results are very similar to the results for the kx scheme shown in Fig. la. 
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Figure 1: a) Dependence of the two-jet exclusive cross section in the kx scheme on the scale 
factor £. The dashed curves are for fi 2 R = fj? F = £ C^2, i Py(*)) 2 - Choosing (Jj. fcy(*)) 2 
as the basic scale yields the solid curves. Choosing Q 2 as the basic scale yields the dotted 
curves. Results are shown for the LO (lower curves) and NLO calculations; b) the same as 
a) for the cone scheme. The kinematical cuts are explained in the text. 

3 The measurement of the gluon density 

Dijet production in DIS at HERA allows for a direct measurement of the gluon 
density in the proton (via "fg — > qq). Some investigations of the feasibility of 
the parton density determination from dijet production hava-been discussed 
in Ref. [6] (see also Ref. [5]). Repond showed first results of the gluon 
distribution determination in NLO in the cone scheme. 

4 Forward Jet Production in the Low x Regime 

Deep inelastic scattering with a measured forward jet with relatively large mo- 
mentum fraction Xj et (in the proton direction) and p^ lab (j) ~ Q 2 is expected 
to provide sensitive information about the BFKL dynamics at low rrEjO. In 
this region there is not much phase space for DGLAP evolution with transverse 
momentum ordering, whereas large effects are expected for BFKL evolution in 
x. In particular, BFKL evolution, is, expected to substantially enhance cross 
sections in the region x « Xj e t LL3'Ea. In order to extract information on the 
ln(l/ir) BFKL evolution, one needs to show that cross section results based on 
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Table 1: Cross sections for n-jet exclusive events in DIS at HERA. See text for details. 



fixed order QCD with DGLAP evolution are not sufficient to describe the data. 
Clearly, next-to-leading order QCD corrections to the DGLAP predictions are 
needed to make this comparison between experiment and theory. 

In Table |l| we show numerical results for the multi jet cross sections with 
(or without) a forward jet. The LD (NLO) results are based on the LO (NLO) 
parton distributions from GRVtHl together with the one- loop (two-loop) for- 
mula for the strong coupling constant. Kinematical cuts are imposed to closely 
model the HI event selection. More specifically, we require Q 2 > 8 GeV 2 , 
x < 0.004, 0.1 < y < 1, an energy cut of E(e') > 11 GeV on the scattered 
electron, and a cut on the pseudo-rapidity rj — — lntan(#/2) of the scattered 
lepton of -2.868 < r){e') < -1.735 (corresponding to 160° < 6(1') < 173.5°). 
Jets arc defined in the cone scheme (in the laboratory frame) with AR = 1 
and \rj(j)\ < 3.5. We require a forward jet with xj et — p z (j)/Ep > 0.05, 
E(j) > 25 GeV, 0.5 < Px(j)/Q 2 < 4, and a cut on the pseudo-rapidity of 
1.735 < r)(j) < 2.9 (corresponding to 6.3° < 8(j) < 20°). In addition all jets 
must have transverse momenta of at least 4 GeV in the lab frame and 2 GeV 
in the Breit frame. 

The cross sections of Table [j] demonstrate first of all that the requirement 
of a forward jet with large longitudinal momentum fraction (xj e t > 0.05) and 
restricted transverse momentum (0.5 < Pr(j)/Q 2 < 4) severely restricts the 
available phase space, in particular for low jet multiplicities. The 1-jet exclusive 
cross section vanishes at LO, due to the contradicting x < 0.004 and Xj e t > 0.05 
requirements. For x << Xj e t, a high invariant mass hadronic system must be 
produced by the photon-parton collision and this condition translates into 



where vtlt and y are the transverse mass and rapidity of the partonic recoil 
system, respectively. Thus a recoil system with substantial transverse momen- 
tum and/or invariant mass must be produced and this condition favors recoil 



with without 
forward jet forward jet 



1 jet (LO) 

2 jet (LO) 

2 jet (NLO) 

3 jet (LO) 



pb 9026 pb 

19.3 pb 2219 pb 

68 pb 2604 pb 

30.1 pb 450 pb 




(1) 
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systems composed out of at least two additional energetic partons. 

As a result one finds very large fixed order perturbative QCD corrections 
(compare 2 jet LO and NLO results with a forward jet in Table |l|). In addition, 
the LO (O(al)) 3-jet cross section is larger than the LO (0(a s )) 2-jet cross 
section. Thus, the forward jet cross sections in Table |l| are dominated by the 
(0(a 2 s )) matrix elements. The effects of BFKL evolution must be seen and 
isolated on top of these fixed order QCD effects. We will analyze these effects 
in a subsequent publication. 
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